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SEVERAL  PROBLEMS  IN  DEVELOPING  A  SATELLITE -BORNE  INFRARED  REMOTE  SENSOR  AND 
WAYS  TO  SOLVE  THEM* 

Gong  Huixing 

Shanghai  Institute  of  Technical  Physics,  Chinese  Academy  of  Sciences 
Received  16  January  1982 

Several  specially-considered  problems  in  developing  a 
satellite-borne  infrared  remote  sensor  are  discussed.  These 
problems  include:  the  definition  of  scanning  period,  the 
scanner  and  its  lubrication  in  outer- space  environment  of 
very  high  vacuum,  and  infrared  radiation  calibration  and  so 
on.  The  ways  to  solve  them  are  suggested. 

In  the  recent  decade,  with  the  appearance  of  a  long-service-life  (free 

from  power  source)  radiation  cooler  and  tellurium-cadmium-mercury  (TeCdHg) 

infrared  detector  with  77RvlOSK  operating  temperature,  detection  rate  D*^ 

10 

1x10  cm. He  /W,  response  time  t<l  is,  and  response  wavelength  X  %  14  trn, 
infrared  remote  sensors  of  optical  mechanical  scanning  type  were  installed 
in  a  long-service-life  operation  meteorological  satellite  and  land  resources 
satellite  to  conduct  periodic  high  discriminabi 1 ity  observations  of  continents, 
atmosphere  and  oceans  of  the  earth.  Since  the  optical  scanning  multiband 
infrared  remote  sensor  has  the  observation  capability  of  the  long  wave  infrared 
wide  band,  the  sensor  occupies  an  important  position  in  space  remote  sensing. 

•All  comrades  in  Section  *720  and  *730  took  part  in  the  work. 
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Figure  1  shows  a  simplified  light  channel  of  a  type  of  multiband 
scanner,  including  a  scan  reflective  mirror,  forming  a  45°  angle  with  the  rotat¬ 
ing  shaft,  which  is  in  the  same  direction  as  satellite  flight.  When  the 
scan  mirror  rotates,  the  instrument  scans  the  direction  perpendicular  to  the 
satellite  flight,  and  receives  radiation  emission  or  reflection  of  the  ground- 
atmosphere  system  by  using  a  fixed  optical  instantaneous  visual  field.  By  taking 
advantage  of  the  forward  motion  of  the  satellite’s  revolution  around  the  earth, 
scenes  of  the  atmosphere,  continents  and  oceans  can  be  obtained  (Fig.  2). 


Fig.  1.  Schematic  diagram  of  light  channel  of  TIROS-N 
AVHRR  four-band  infrared  scanner. 

Key:  (a)  Motor;  (b)  Target  radiation;  (c)  Dispersion 
lens;  (d)  Detector. 


This  paper  gives  a  brief  discussion  concerning  problems  and  solution 
methods  in  the  development  of  this  type  of  remote  sensor. 

1.  Determination  of  scanning  speed 

In  order  to  maintain  a  constant  scale  (with  respect  to  time  and  location) 
of  the  received  scene  photographs,  the  scanner  requires  that  the  satellite 
orbit  be  a  circle.  In  that  case,  the  moving  speed  of  the  satellite  projection 
on  the  ground  surface  is: 
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irnrx  *+*  • 


In  the  equation,  R  is  the  earth's  radius;  h  is  the  orbit  altitude;  and  the 
Kepler  constant  u *3.986x10^  km^/s 


Fig.  2.  Diagram  showing  the  principle 
of  photographing  an  earth  scene  by  a 
scanner. 

Key:  (1)  Satellite;  (2)  Orbit;  (3)  Earth; 

(4)  Scanning  lines;  (5)  Satellite  projection 
locus;  (6)  Direction  of  earth's  rotation. 


If  the  latitude  of  the  instantaneous  visual  field  of  a  detector  (of  the 
instrument)  is  a  along  the  satellite  flight  direction,  in  order  to  avoid 
misscanning  of  scenes  on  the  ground  surface,  the  scanning  period  T  of  the 
scanner  should  be: 


T 


N  'a'h(R+h) 


;  r+\ 

V  n  • 


(2) 


In  the  equation,  N  is  the  number  of  elements  of  the  detector. 


2.  Requirements  on  orbit 


In  the  radiation  sounding  of  the  infrared  channel  (iO.S~12.S,*.m)  of  an 
atmospheric  spectral  window,  a  TeCdMg  detector  is  required;  the  detector  is 
installed  on  a  second-stage  cold  block  of  the  two-stage  radiation  cooler  with 
an  operating  temperature  of  10SK.  In  order  to  maintain  sufficient  exchange  of 
heat  radiation  between  the  radiation  cooler  and  the  3K  cold  space  environment, 
the  opening  of  the  radiator  should  aim  consistently  at  the  cold  space. 

In  a  satellite  with  three-axis  attitude  stability  with  respect  to  the 
ground,  the  drift  and  rolling  shafts  of  the  satellite  rotate  around  the  pitch 
axis  at  the  same  angular  velocity  as  the  satellite's  revolution  around  the 
earth.  Since  the  pitch  axis  is  perpendicular  to  the  orbital  plane,  the  axis  is 
relatively  stable  in  space.  Hence,  in  the  polar  region  and  the  synchronizing 
orbit  of  the  earth,  the  opening  of  the  radiation  cooler  is  toward  the  shady 
direction,  the  same  as  the  pitch  axis.  In  order  to  avoid  sunlight  shining 
onto  the  opening  of  the  radiation  cooler,  it  is  necessary  to  maintain  generally 
a  constant  angle  between  the  sunlight  and  t-he  orbital  plane  of  the  satellite; 
this  is  the  solar  synchronizing  orbit.  The  method  of  synchronizing  with  the 
sun  is  to  control  the  launch  dip  angle  i  of  the  orbit.  Mien  i  satisfies 

i-oo.->[-0.0985(l-e»)*(«)W]  (8) 

(e  is  the  orbit  eccentricity;  and  a  is  the  semi-major  axis  of  the  orbit),  the 
perturbation  of  the  satellite's  motion  by  the  protruding  mass  of  the  equator 
can  allow  the  rotation  of  the  orbital  plane  around  the  earth's  rotational 
axis  to  cancel  the  variation  between  the  sunlight  and  the  orbital  plane 
because  of  the  earth's  revolution  around  the  sun.  In  tne  synchronizing  orbit 
around  the  sun,  the  annual  variation  of  the  intersecting  angle  between  the  sun 
and  the  orbital  plane  does  not  exceed  +23.5°. 

5.  Drive  device  of  scanning  mirror 

The  earth-scene  signals  used  for  satellite  photography  can  be  converted 
into  pictures  through  ground  receiving  stations  by  machine  synchronization. 


When  satellite  picture  transmission  uses  a  single  channel  simulated  signal  for 
output,  the  longterm  and  instantaneous  stability  of  the  rotational  speed  of 
the  scanning  mirror  should  be  strictly  controlled  to  avoid  slanting  and  fluctua¬ 
tion  of  pictures  received  by  the  ground  because  the  synchronized  reference 
signals  of  rotational  speed  (of  the  scanning  mirror)  do  not  accompany  the  output 
signal. 

In  order  to  simplify  the  transmission  structure,  to  solve  the  wearing  and 
lubrication  problems  of  the  gears,  and  to  increase  service  life  and  reliability, 
a  TDB-500  low-speed  permanent -inagnetic-type  ac  synchronous  motor  is  used  to 
directly  drive  the  45°  scanning  mirror.  The  TDB-500  is  the  axial -direct ion 
exciting  magnet,  operating  on  the  principle  of  varying  gas  gap  magnetic  conduct¬ 
ance  due  to  open  slots  of  the  stator  and  rotor.  The  rotational  speed  N  of  the 
rotor  is: 

N=60x  (frequency  of  power  source/number  of  open  slots  of  the  rotor) 
(rpm)  .  (4) 

The  number  of  open  slots  of  the  rotor  is  50;  when  the  frequency  of  the  power 
source  is  100  Hz,  the  motor  rotational  speed  is  120  rpm.  Since  this  is  a 
synchronous  motor,  the  average  rotational  speed  of  the  rotor  is  determined  by 
the  frequency  of  the  power  source.  By  using  a  power  source  with  quartz 
oscillation  frequency  division,  a  frequency  stability  of  3xl0*6  can  be  attained. 
The  motor  operates  on  the  speed  reduction  principle  of  an  electromagnet  with 
50  times  the  specific  synchronizing  moment  of  an  ordinary  motor.  Through 
measurement,  the  fluctuating  angle  of  the  scanning  mirror  is  less  than  0.8  mrad. 

The  greatest  flaw  in  using  a  synchronous  motor  is  that  it  is  relatively 
complex  in  starting  up  because  the  greatest  synchronizing  angle  is  i.8°  for 
a  motor  with  a  speed  reduction  ratio  of  50,  When  the  rotational  speed  is  100 
rpm,  it  is  equivalent  to  starting  up  the  rotor  component  with  an  inertia  (J= 

*  *y 

i.8xl0  g-cnr)  within  a  duration  of  5  ms.  Based  on  the  relationship  of  the 
•  ^  2 

starting  moment  M*2r “J (N/60)  x  speed  reduction  ratio,  the  motor  should  have  a 
moment  of  T^xlO5  g-cm,  which  considerably  exceeds  the  moment  output  capacity 
of  the  TDB-500  motor.  Hence,  a  variable  frequency  scheme  is  selected.  Under 
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the  situation  of  constant  motcr  voltage  of  the  power  source,  when  the  source 
frequency  is  reduced,  the  stator  current  increases  and  the  starting  moment 
rises;  the  inertia  moment  of  the  rotor  decreases  with  the  square  law.  Figure  3 
shows  that  the  motor's  starting  moment  is  equal  to  its  inertia  moment  when 
the  rotational  speed  is  approximately  25  rpm.  The  frictional  moment  of  the 
rotor  component  is  ISOg-cm;  the  maximum  synchronizing  moment  is  900g-cm  when 
the  motor  has  a  source  voltage  of  45V;  and  the  initial  rotational  speed  is  10  rpm 
for  variable-frequency  starting.  From  Fig.  3,  the  allowance  of  the  motor 
starting  moment  is  9,  and  the  allowance  of  the  operating  moment  is  5;  this 
ensures  operational  reliability  of  the  motor.  In  order  to  solve  the  problem 
of  wide-range  f .  quency-varying  speed-adjustment  motor  starting,  a  numerical 
phase  splitting  variable  frequency  power  source  was  developed. 


/JlPfr*  rm) 


Fig.  3,  Startim,  characteristics  of  scanning  mirror 
component . 

Key:  (a)  Moment;  (b)  Inertia  moment;  (c)  Starting 
moment;  (d)  Rotational  speed. 


4.  Lubrication  of  rotating  components  in  very  high  vacuum 


Above  a  700-km  orbit,  the  degree  of  vacuum  of  a  satellite’s  space 

-8 

environment  is  higher  than  10  Torr.  For  an  optical  scanning  system,  there 
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are  problems  of  lubricating  rotational  components  and  of  operational  life 
span  in  conditions  of  very  high  vacuum. 


There  are  two  generally  used  space  lubrication  methods:  solid  lubrication 
and  oil-gas  lubrication.  According  to  availability  abroad  and  existing 
domestic  conditions,  we  selected  an  oil-gas  lubricating  structure  with  oil 
bearing  retainer  frame  and  labyrinth  sealing  of  the  protruding  part  of  the 
rotational  shaft.  This  plan  was  tested  for  long  service  life  under  high- 
vacuum  ground  conditions. 

Analyzing  the  ring-shaped  pipeline  in  Fig.  4,  r^  represents  the  inner 
diameter  of  the  rotational  shaft;  r-  represents  the  outer  diameter  of  the  shaft 
hole.  When  the  average  free  stroke  X  of  gas  molecules  is  considerably  greater 
than  the  shaft  radius,  the  gas  conduction  F  of  the  ring-shaped  pipeline  is: 

F- 80.48^5  <*?-<< '•—'*!  .  [L,s]  (5) 

In  the  equation,  T  is  the  absolute  temperature;  M  is  the  gram  molecular  weight 
of  the  gas;  and  j,  is  the  length  of  pipeline. 


Fig.  4.  Ring-shaped  pipeline. 

The  gas  flow  U  fin  the  ring-shaped  pipeline1)  is: 

r  :/•,  -r:.F,  -Torr-I.  *],  (6> 

In  the  equation,  fP^-P,)  is  the  pressure  difference  between  two  sides  of  the 
ring-shaped  pipeline. 


From  Eq.  (6),  the  gas  outflow  from  the  ring-shaped  pipeline  is  related  to 
pressure  difference  and  gas  conduction.  In  outer  space,  P^O,  therefore  the 
pressure  difference  between  two  sides  of  the  ring-shaped  pipeline  is  primarily 
determined  by  the  gas  pressure  within  the  rotational  device.  Based  on  the 
analysis  mentioned  above,  only  by  selecting  a  lubricating  oil  with  low 
saturated  vapor  pressure,  and  by  controlling  the  length  of  the  gap  of  the 
protruding  portion  of  the  rotational  shaft,  and  by  storing  sufficient 
lubricating  oil  in  the  motor  oil  storage  device  and  the  bearing  retainer  frame 
can  a  saturated  vapor  pressure  of  the  lubricating  oil  be  maintained  for  a  long 
time  in  the  motor  to  achieve  the  purpose  of  bearing  lubrication. 

We  selected  #114  silicone  oil  as  lubricant  with  a  gram  molecular  weight  of 
8047,  lxlO7  Torr  as  the  saturated  vapor  pressure  at  50°C,  and  dimensions  of 
labyrinth  clearance:  ijsl.l  cm,  r^l.095  cm,  and1=l  cm.  As  calculated  from 
Eq.  (6),  the  total  amount  of  lubricating  oil  (escaping  into  outer  space  from 
two  openings  of  labyrinth  clearance  in  a  year)  is  1  mg.  If  100  mg  of  lubricating 
oil  is  stored  in  the  scanner,  two  years’  continuous  operation  of  the  scanner 
in  orbit  can  be  ensured. 

In  order  to  r  rtify  the  above-mentioned  scheme,  we  conducted  a  simulation 

test  by  placing  th  •  Tanning  components  into  a  KM-0.5  container  and  operated  a 
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Model  2TL-500  spattering  titanium  pump  up  to  3^6x10  Torr.  Tne  scanner  operated 
continuously  for  half  a  year  at  120  rpm;  no  appreciable  variation  of  bearing 
friction  moment  was  discovered. 

5.  Radiation  calibration  in  flight 

In  order  to  deduce  radiation  intensity  and  temperature  of  a  target,  from 
the  output  signal  of  an  infrared  scanner,  it  is  required  to  conduct  radiation 
calibration  of  the  scanner. 

With  a  cold  space  background  (3K)  in  outer  space,  the  radiation  emissivity 

I  *)  C 

is  only  1/1C  4"  of  a  normal-temperature  target;  this  can  be  the  :ero  radiation 
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basic  reference  standard  of  the  scanner.  The  amplitude  of  the  oucput  signal  of 

an  infrared  TdCdHg  detector  is  of  microvolt  magnitude.  In  order  to  avoid  dc 

drift  in  the  amplification  process  of  weak  signals,  ac  ccupling  is  adopted  for 

the  infrared  prime  amplifier.  Since  the  coupling  capacitor  separates  the  dc 

component  into  signals,  the  output  of  the  prime  amplifier  only  reveals  the 

variation  in  radiation  received.  Ir.  order  to  measure  the  absolute  intensity, 

when  a  scanner  observes  cold  space,  a  logic  control  circuit  produces  a  dc 

restoration  pulse  to  operate  a  dc  restoration  capacitor  C  of  a  restoration 

amplifier  to  rapidly  discharge  against  the  ground  via  a  conduction  diode  D, 

thus  establishing  a  zero  signal  power  level  (see  Fig.  5)  corresponding  to  zero 

radiation  in  space.  The  diode  is  shut  off  after  the  dc  restoration  pulses 

pass.  If  RC>(l/Ul)  (w  is  the  lower  frequency  limit  of  the  channel  amplifier), 
l  L 

the  capacitor  C  exerts  a  clamp  function  of  zero  power  level.  Later,  the 
amplitude  of  the  output  signal  is  proportional  to  the  absolute  radiation 
intensity  of  the  target  received  by  the  scanner. 


Vx-  —KQA r0(X)P (X,  Tm)dX, 
*  (*) 
P(x* 


Key:  *)  target. 

-2  2 

In  the  equation,  ^target  is  the  target  temperature;  0^=3.74x10  W-cm  ;  Cy 
1.438  cnrK;  tq(A)  is  the  optical  efficiency;  R(A)  is  the  response  rate  of  the 
detector;  A  is  the  radiation  receiving  area  of  the  optical  system;  ft  is  the 
instantaneous  visual  field;  and  K  is  the  electronic  amplification  time. 


In  order  to  measure  the  radiation  absolute  intensity  of  a  target  correspond 
ing  to  signal  V,  a  reference  black  body  (for  example,  the  shell  of  an  instrument 
can  be  made  into  a  reference  black  body)  is  introduced  in  an  instrument.  If  the 
output  signal  is  when  the  scanner  observes  the  reference  black  body,  then: 

i  •».  (*) 

^(k)T,(X)P(/.,  T.idX,  (8) 


Key:  *)  reference. 


Fig.  5.  Diagram  showing  principle  of  dc 
restoration  amplifier. 

Key:  (1)  TeCdHg;  (2)  Prime  amplifier; 
(3)  Amplification;  (4)  Output;  (5)  dc 
restoration  pulse. 


In  the  equation,  Treference  *s  the  temperature  of  a  reference  black  body,  then 


O) 

v  l  :#<X)t,A)P<X  T.Vtt. 


<») 


Key:  (1)  target;  (2)  reference. 

In  the  above  equation,  the  ratio  of  and  can  be  measured  at  the  ground 
surface,  ^reference  can  he  obtained  by  temperature  remote  sensing  of  the 
reference  black  body.  Hence,  the  equivalent  black-body  temperature  of  the 
target  can  be  revealed  by  using  Eq.  (9). 


6.  Radiation  calibration  of  ground  surface 


The  response  rate  R(x)  of  a  detector  can  be  kept  unchanged  by  accurately 
controlling  cold-block  temperature  of  the  radiation  cooler.  However,  the 
optical  efficiency  tq(X)  of  the  scanner  will  vary  with  the  instrument 
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temperature  (primarily  determined  by  the  color  division  lens  and  light 
filtering  lens).  Limited  by  the  structure  and  weight  of  the  satellite-borne 
instrur'  the  blackness  of  the  reference  black  body  cannot  be  made  very 
high.  Hence,  radiation  calibration  should  be  done  on  the  ground  before 
launching  of  the  scanner  in  order  to  obtain  a  set  of  the  following  calibration 
curve  for  users'  reference  (see  Fig.  6). 


B  *3  T. 

Fig.  6.  Radiation  calibration  curves  of 
radiometer. 

Key:  (1)  Direction  of  temperature  increase 

of  the  shell;  (2)  Equivalent  black  body 
temperautre  of  the  target. 

For  a  scanner  including  a  radiation  cooler,  the  radiation  calibration  of 
the  infrared  channel  should  be  conducted  in  a  vacuum  system  of  simulated  cold 
background.  The  temperautre  of  the  cold  background  should  be  below  33K;  in 
that  case,  the  effect  on  the  cooling  capacity  of  the  second-stage  cold  block 
is  less  than  1  percent.  During  the  experiment,  the  radiation  cooler  is  aimed 
at  the  cold  shield,  and  a  liquid  nitrogen  heat  sink  (generally,  the  temperature 
is  90K;  the  radiation  intensity  within  10.5a,12.5  um  is  only  1/10,000  of  the  270K 
black  body  target)  of  a  vacuum  container  is  used  as  the  reference  datum  of  zero 
radiation  during  dc  restoration  of  signals. 

During  calibration,  a  honeycomb  surface  source  can  be  used  as  the  standard 
black  body;  dimensions  of  the  surface  source  should  he  able  to  fill  the  light 
passage  diameter  and  instantaneous  visual  field  of  the  scanner.  Based  on 
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manuscript  reports,  if  for  each  honeycomb  element  of  6  mm  wide,  25  mm  long,  and 
0.05  mm  in  wall  thickness,  coating  black  paint  on  the  surface  with  thermal 
radiation  efficiency  £=0.95,  then  the  equivalent  thermal  radiation  efficiency 
Cj=0.996.  This  corresponds  to  a  temperature  difference  of  0.23K  (at  270K) 
with  an  ideal  black  body. 


In  order  to  further  increase  the  thermal  radiation  efficiency  of  the  surface 
source  black  body,  the  honeycomb  surface  source  board  can  be  made  into  a  cavity, 
as  shown  in  Fig.  7.  The  ratio  between  diameter  and  length  of  the  opening  should 
be  1:2.  According  to  Gouffe's  formula,  we  estimate  the  equivalent  thermal  radia¬ 
tion  efficiency  e,  of  the  cavity  body  as: 


g,(l  -t-m1) 


.Si-  .4  S‘ 


(10) 


Fig.  7.  Cavity-shaped  honeycomb 
surface  source  black  body. 


Ill  the  equation,  at  the  cavity  wall  the  equivalent  thermal  radiation 
efficiency  (of  the  cavity  wall)  Cj=0.996;  the  ratio  (between  opening  area  and 
the  cavity  surface  area)  A/S=l/12, 

»-(l- '»)(-*  6.)-0.0Ul(-i.--L)-6:<tU  \ 

is  the  surface  a.ea  of  a  sphere  with  its  diameter  equal  to  the  cavity  depth. 
Then  we  can  calculate:  C2=0--9^"6-  Th*5  reveals  that  the  error  of  equivalent 
black  body  temperautre  (at  270K)  is  less  than  0.01K  due  to  thermal  radiation 
efficiency.  Then  the  radiation  calibration  accuracy  of  the  scanner  is  determined 
by  the  temperature  error  of  the  cavity-shaped  surface  source  black  body,  and 
the  signal-to-noise  ratio  and  signal  processing  method  of  the  output  signals 
of  the  scanner. 
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The  effoots  of  the  surface  parameters,  <  end  T„  on  atmospheric  temperature  retrie¬ 
vals,  end  the  possibility  of  reoovering  simultaneously  both  c,  T,  and  atmospherio 
temperature  profiles  are  studied.  The  rosults  of  numerioal  simulations  show  that,  when 
one  of  the  two  surface  parameters  Is  known,  atmospherio  temperature  profiles  oan  be 
derived  end  et  the  same  time  the  other  surface  parameter  can  be  determined  accurately; 
and  even  If  both  a  and  T,  are  unknown,  atmospherio  temperature  profiles  and  the  two 
parameters  oan  still  be  recovered  with  good  accuracies. 


I.  FOREWORD 

Satellite  remote  sensing  can  provide  the  atmospheric  temperature  contour 
lines  from  the  ground  surface  to  the  top  of  the  stratosphere  over  the  entire 
globe  (including  continents  and  oceans).  Since  the  ground  characteristics 
are  more  complicated  than  those  of  oceans,  the  thermal  radiation  efficiency 
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and  ground  temperature  are  not  easy  to  predict;  correspondingly,  retrievals 
of  the  atmospheric  temperature  over  land  are  difficult. 

For  a  group  of  temperature  measuring  channels  of  a  satellite,  the  measured 
upward  radiation  intensity  of  the  i-th  channel  can  be  expressed  as: 

R,- [€tB,(.T.)  +  +  \  B,(T)dTt(P),  (1) 

(«-t,  a,  M) 

Here,  M  represents  the  number  of  channels;  e  is  the  thermal  radiation  efficiency 
at  the  ground  surface;  B  is  the  Planck  function;  R|  is  the  effective  downward 
atmospheric  radiation  intensity;  t(P)  is  the  transmissibility  from  the  atmo¬ 
spheric  pressure  P  to  the  top  of  the  atmosphere;  the  subscript  i  represents 
the  amount  of  the  i-th  channel;  T  is  the  atmospheric  temperature  at  the 
atmospheric  pressure  P;  and  Ts  and  Ps  represent,  respectively,  the  temperature 
and  atmospheric  pressure  at  the  ground  surface. 

From  the  above  transmission  equation  (1),  to  solve  for  the  atmospheric 

temperature  contour  lines,  it  is  assumed  that  the  efficiency  of  the  ground 

thermal  radiation  is  known  (generally,  the  ground  surface  is  considered  as  a 

black  body;  that  is,  c*l),  assuming  that  Ps  is  the  average  atmospheric  pressure 

at  the  sea  surface  and  T  is  measured  from  the  channel  of  the  spectral  window; 

s  f  1  31 

all  these  parameters  are  known1  J.  Or  T  is  considered  unknown  and  is 

s  f 4  si 

derived  together  with  the  atmospheric  temperature1  ’  J.  For  the  sea  surface, 
the  thermal  radiation  efficiency  is  close  ti  1,  since  the  thermal  radiation 
efficiency  can  be  considered  as  a  known  constant  at  least  for  a  selected  group 
of  remote  sensing  channels.  Hence,  the  above  processing  is  appropriate.  There 
are  different  conditions  over  the  land;  although  the  atmospheric  pressure  Ps 
at  the  ground  surface  can  be  determined  based  on  the  topographic  altitude  of 
various  sites,  or  directly  measured^  by  laser  radar  from  a  satellite,  yet  the 
thermal  radiation  efficiency  varies  obviously  with  different  characteristics 
of  the  ground  surface,  and  numerically  the  value  of  efficiency  nay  be  consid¬ 
erably  different  from  1.  Moreover,  the  value  of  thermal  radiation  efficiency 

;  is 
t 


is  different  for  different  wavelength  values;  obviously,  the  value  is  not  a 

constant:  Since  the  thermal  radiation  efficiency  at  the  ground  surface  is 

unable  to  be  accurately  determined,  the  value  also  affects  the  ground  surface 
fllMSl 

temperature1  measured  by  use  of  channels  at  the  spectral  window;  thus, 

the  derived  Ts  may  have  errors.  Hence,  the  retrievals  of  the  atmospheric 
temperature  over  land  are  more  difficult  than  for  the  sea  surface.  It  is  more 
sc  for  the  non-statistical  retrieval  method. 

At  present,  the  N  series  TIROS  [television  and  infrared  observation 
satellite]  satellites  of  the  United  States  apply  the  regression  method  of  zonal 
statistics  for  temperature  retrievals.  In  this  method,  the  entire  globe  is 
divided  into  five  latitude  zones:  one  zone  from  30*S  to  30*N,  and  two  zones 
each  from  the  pole  to  60*,  and  from  60*  to  30*.  Then  the  respective  regression 
coefficients  can  be  derived.  These  regression  coefficients  can  be  continuously 
corrected^1^  by  using  consistent  data  values  from  radio- sondages  and  satel¬ 
lites.  On  the  ether  hand,  researchers  are  studying  ways  of  using  other  than 
the  conventional  data;  the  accuracies  of  thesd  kinds  of  retrieval  methods  are 
comparable  to  that  using  the  conventional  data. 

Authors  of  the  manuscript  [17]  studied  the  retrieval  method  described  in 
manuscript  [4],  considering  the  non-black-body  ground  surface  as  a  black  body, 
and  studying  the  retrieval  effects  of  temperature  and  atmospheric  pressure 
at  the  ground  surface.  In  numerical  experiments,  the  given  values  of  the 
thermal  radiation  efficiency  at  the  ground  surface  are  0.85  and  0.95  at  the 
4.3  urn  band,  and  0.95  and  0.98  at  the  IS  ym  band.  As  revealed  by  results,  the 
errors  cannot  be  neglected. 

The  paper  studies  the  effects  of  c  and  Tg  on  temperature  retrievals, 
and  investigates  the  possibility  of  simultaneously  deriving e  and  Tg  on  one 
hand,  and  the  atmospheric  temperature  on  the  other. 

II.  MATHEMATICAL  PROCESSING 

The  15  urn  CO,  band  is  selected  as  the  temperature  measuring  channel; 
then  on  the  right  side  of  equation  (1),  the  ground  surface  reflection  term  of 
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space  radiation  is  nuch  smaller  than  the  ground-surface  emission  term. 
Therefore,  no  appreciable  errors  will  be  introduced  by  using  appropriate 
modes  for  approximation.  As  derived  from  manuscript  [17],  the  effective 

A 

downward  atmospheric  radiation  intensity  R^|  can  be  expressed  as 

(2) 

In  the  equation,  the  factor  can  be  considered  as  a  constant  related  to  the 
channels;  the  value  of  Fj^  varies  very  little  in  different  temperature  layers 
of  the  atmosphere.  Therefore,  a  line-by-line  calculation  method  can  be  used 
in  advance  to  derive  a  group  of  constants  [F!]  corresponding  to  a  group  of 
remote-sensing  channels.  Thus,  equation  (1)  can  be  rewritten  as: 

JJ.-{*.+P,,(l-£,)[l-T,(P.)]>B,(r.)T,(P.)  +  f1  (3) 

It  is  assumed  that  the  atmospheric  pressure  Ps  at  the  ground  surface  is  known 
in  the  wave  band  for  a  selected  group  of  channels,  the  ground  surface  can  be 
considered  a  gray  body  (i.e.,  Ej*e»  i*l,  2,  ...,  M) ;  then  the  following 
equation  can  be  derived  from  dq.  (3): 

R-RS  +  AR?',  W 

In  the  equation, 

i*;“- [<a>+i:(l-fu>)(l-r?)]Si(Tnr<:,+  P  B,(Tu>)drr(P),  (5) 

AR"  -  A  -  Rtl>  -  [1  -  (1  -  ri?)  ]  5,(770 Ta° & “ 

n  (i  -  *"0  (I  -  t(°)  ]  r'»  1  or? 

Ol ,  I  it* 

<6) 

In  equations  (5)  and  (6),  t,  5t.(PJ  (as  in  the  following).  The  entire  atmo- 
spheric  layer  above  the  ground  surface  can  be  divided  into  N  thin  layers;  by 
using  the  numerical  integration  formula,  the  equation  (6)  can  be  expressed  as 

j r:  -  -  R,n  -  .47 j a',".  <»  - 1, 2,  -m)  co 

In  the  equation. 
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A$-^p-\r„  WPi-i)-r?Wh  0  -1.  2.  N) 
J&M  -  0'°+^  (i-«">)  (1  -  03  <!'  — ^  P- j 


-  n  a  -  o]«.(T?»)ti?# 

J.v;  -JT5«>,  (j-1.  2,  ...N), 

•  -  JTJ\ 

JA\4V.«  Jt"», 

P0aO  and  *r^(Po)sl;  is  the  average  temperature  of  the  j-th  layer  of  the 
atmosphere.  By  using  the  method  in  manuscript  [S] ,  equation  (7)  can  have 
solutions  in  matrix  form: 


JX">  -  (  A‘J1)  *  [A">(  A(,))  *  +  rl]  -MR"-,  (8) 

X«“>-X(',  +  dX,,'(  (B) 

In  equation  (8),  I  is  the  MxM-dimension  unit  matrix;  r  is  a  smoothing  factor; 
the  superscript  represents  times  of  successive  substitution;  and^=0  is  the 
quantity  corresponding  to  the  initial  value.  The  process  of  successive 
substitution  stops  at 

Wl<«,  (»-l,  2,  ...M) 

and  c  is  the  observation  error. 

III.  RESULTS  AND  DISCUSSIONS  ON  NUMERICAL  TESTS 

The  author  and  his  colleagues  used  channels  (see  Table  1)  of  two  groups 
of  triangular  response  function  described  in  manuscript  [5],  The  numerical 
experiments  of  temperature  retrievals  were  conducted  by  assuming  in  the 
following  (different  gro  'nd-surface  conditions)  that  the  atmospheric  pressure 
at  the  ground  surface  is  1000  mb;  the  measurement  error  is  0.3  erg* s"1 -cm  *• 
ster'l;  and  the  thermal  radiation  efficiency  (at  the  ground  surface)  varies 


between  0.8  and  0.99. 
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the  coefficient  matrix 


1.  e  Md  T  are  known;  then  AT,ihr=to, 

fti  s 

Av  degenerates  into  a  MxN-dimer.sion  matrix. 

ft-) 

2.  It  is  assumed  that  the  ground  surface  is  a  black  body  (c  sO),  and 
Ts  is  an  unknown.  Then,  Ae^sO;  the  coefficient  matrix  A^  is  a  matrix 

of  Mx(N+l) -dimension. 

3.  e  is  a- known  quantity  but  T  is  unknown.  Then, Ac  =0;  the  coefficient 

ftl  s 

matrix  A  is  a  matrix  of  Mx(N-*l)  dimensions. 

4.  e  is  an  unknown  quantity  but  Tg  is  known 
deviations  of  0®C,  2®C  and  5*C  between  Tfi  and  the 
ture.  Thos,  At^sO;  the  coefficient  matrix  A^ 

Mx(N+l) -dimension. 

5.  Both  £  and  T  are  considered  unknown  quantities  to  be  determined. 

s  d  1 

Then,  the  coefficient  matrix  Av  ‘  is  a  matrix  of  Mx(N*2) -dimension. 

Table  1.  Number  of  channels  and  center  frequencies  of  various  channels  for  two 
groups  of  temperature  measurement  channels. 


fa)*  9 

_ (e)  * 

a  a  -t-  «  * 

(dJ7 

608,  6:0.3.  692.  704.5,  717,  732,  747 

Cc) 

17 

668.  672,  677.  682.  ( 

•87,  692,  097,  700.3,  702.  707,  712.  717, 

722,  732.  747,  752,  757 

by  assuming  that  there  are 
"real"  ground- surface  tempera- 
degenerates  into  a  matrix  of 


Key:  (a)  Sequence  number;  (b)  First  group;  (c)  Second  group;  (d)  Number  of 
channels;  (e)  Center  frequencies  (cm-*)  of  various  channels. 


Table  2  shows  the  results  obtained  in  numerical  experiments  in  different 
conditions  at  the  ground  surface.  The  results  include  the  root  mean  square 
R>!5(*C)  of  the  atmospheric  temperatures  of  37  pressure  layers  in  all  from 
10  mb  to  the  ground  surface;  the  root  mean  square  deviation  (*C)  from  the 
ground- surface  temperature  of  19  layers  (of  the  37),  and  the  root  mean  square 
deviation  Ae  of  the  thermal  radiation  efficiency  of  the  ground  surface. 
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Remarks:  (1)  The  column  marked  with  *  shows  results  of  210  numerical  experiments;  the  other  columns  show 

the  results  of  2«9  numerical  experiments.  (2)  c  and  Ts  are,  respectively,  the  ground -surface  thermal 
radiation  efficiency  and  the  ground-surface  temperatures;  c’  and  T'  are,  respectively,  the  assumed  known 
ground-surface  thermal  radiation  efficiency  and  the  ground-surface  temperature. 

Key:  (a)  Number  of  layers;  (b)  Pressure;  (c)  Number  of  channels  and  ground-surface  conditions;  (d)  Unknown 
(e)  Ground  surface. 


From  Table  2,  results  obtained  by  use  of  17  channels  are  obviously  better 
than  those  obtained  for  7  channels;  this  is  consistent  with  results^  of 
retrievals  of  the  atmospheric  temperature  over  the  sea  surface.  Let  us 
investigate,  in  the  following,  the  results  obtained  by  use  of  17  channels  of 
temperature  measurement.  Under  the  condition  that  both  ground-surface  param¬ 
eters,  e  and  Ts,  are  accurately  known,  the  retrieval  results  of  the  atmospheric 
temperature  are  the  best.  However,  by  only  having  one  of  two  (t  and  T  )  as  a 
known  quantity,  we  also  obtained  the  same  good  results;  moreover,  at  the  sane 
time,  another  ground-surface  parameter  can  be  accurately  obtained.  When  Ts 
is  known,  the  obtained  (root  mean  square)  deviation  of  the  thermal  radiation 
efficiency  is  0.0022;  when  e  is  known,  the  obtained  (root  mean  square)  deviation 
of  Tg  is  0.1S#C.  However,  when  errors  develop  in  one  ground-surface  parameter, 
other  errors  also  exist  in  the  retrieval  of  the  other  ground- surface  parameter 
(and  at  the  same  time,  the  error  is  increased  in  retrieval  of  the  atmospheric 
temperature).  This  effect  exists,  primarily,  in  several  layers  close  to  the 
ground  surface. 

Also  from  Table  2,  when  Tg  and  e  are  known,  the  retrieval  accuracy  of  the 
atmospheric  temperature  is  also  comparable  to  the  condition  that  both  Ts  and 
e  are  known.  In  addition,  values  of  Tg  and  e  can  be  simultaneously  and 
accurately  obtained;  the  deviations  of  the  root  mean  square  of  the  two  are, 
respectively,  1.13°C  and  0.155.  While  one  mode  of  atmosphere  and  two  radio- 
sondage  temperature  contour  lines  are  considered  the  "real"  situation.  Figure  1 
shows  the  comparison  of  retrieval  results  assuming  that  both  T  and  e  are  known, 
and  both  Tg  and  e  are  unknown  in  two  different  conditions  at  the  ground  surface. 
Of  the  retrieval  results  under  the  two  conditions,  there  are  only  small  differ¬ 
ences  at  layers  close  to  the  ground  surface;  these  two  results  are  almost 
entirely  identical  for  higher  layers. 

While  accurately  knowing  the  contribution  of  the  ground  surface  to  the 
emergent  radiation,  there  is  no  effect  on  retrieval  of  the  atmospheric 
temperature  by  the  ground-surface  parameters.  The  situation  is  similar  when 
both  the  ground-surface  parameters  t  and  Tg  are  known.  Under  other  situations, 
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there  is  some  difference  from  the  above-mentioned  situation  for  the  retrieval 
results  of  the  atmospheric  temperature;  this  difference  is  primarily  due  to 
errors  in  ground-surface  parameters.  However,  relatively  speaking,  errors  in 
ground-surface  parameters  do  not  have  an  appreciable  effect  on  differences 
in  the  atmospheric  temperature.  From  Table  2,  the  effect  of  errors  (between 
e  and  T  )  is  close  to  a  proportional  relationship.  In  Equation  (5),  the  main 
contribution  to  the  ground  surfu.ce  is  from  the  ground-surface  emission  term; 
this  also  nea.'r.  "hat  the  functions  of  e  and  Tg  are  mutually  compensated  to  a 
very  considerable  extent.  This  indicates  that  there  is  a  closer  relationship 
between  two  ground-surface  parameters  than  between  two  _  rameters  and  the 
atmospheric  temperature.  Just  because  both  inherent  relation  and  relative 
independence  exist  between  the  ground-surface  parameters  ar.d  the  atmospheric 
temperatures,  we  can  possibly  determine  ground-surface  parameters  with 
satisfactory  accuracy  while  deriv  lag  the  atmospheric  temperatures. 

We  know  that  following  the  correction  of  the  atmospheric  effect,  when  one 
of  the  two  (ground-surface  parameters,  e  and  T  ,  of  a  channel  of  the  spectral 
window)  is  known,  the  other  parameter  can  be  determined.  Both  r-  and  Ts  can 
be  simultaneously  determined  for  two  channels  of  the  spectral  window  with  the 
same  thermal  radiation  efficiency.  Under  situations  of  this  type,  errors  of 
the  obtained  ground- surface  parameters  are  determined  by  whether  or  not 
corrections  of  the  atmospheric  effect  are  accurate  if  other  errors  are  not 
considered.  \nong  temperature-measurement  channels  we  used,  those  channels 
having  the  ground-surface  transmissibility  not  equal  to  zero  also  include 
signals  of  ground-surface  parameters.  However,  these  channels  are  within  the 
absorption  zone,  with  considerable  effect  played  by  the  atmosphere,  so  it  is 
even  more  difficult  in  correction.  Comparing  the  obtained  temperature 
contour  lines  with  the  real  contour  lines  in  temperature  retrievals,  oftentimes 
when  positive  deviations  exist  in  some  layers,  negative  deviations  exist  in 
the  adjacent  layers.  Therefore,  the  contribution  of  the  entire  atmosphere 
(the  last  term  on  the  right  side  of  equation  (3))  is  still  possibly  close  to 
the  real  value.  For  the  ground- surface  transmissibility,  there  is  no  major 
effect  on  different  temperature  layers  of  the  atmosphere.  Thus,  equations  of 
ground-surface  parameters  can  be  derived  from  equation  (3);  it  is  possible 
to  determine  c  and  T$  like  that  in  the  spectral  window  zone.  The  number  of 
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signals  of  ground-surface  parameters  contained  by  these  channels  is  smaller 
than  those  in  the  spectral  window  zone;  errors  with  corrections  of  the 
atmospheric  effect  may  be  relatively  great.  Furthermore,  by  considering  the 
existence  of  the  observation  errors,  the  ground-surface  parameters  cannot  be 
determined  as  well  as  that  for  the  situation  of  the  spectral  window  zone 
with  the  same  number  of  channels.  For  compensation  of  this  point,  it 
is  re^uir.-d  to  adopt  more  channels  including  signals  of  the  ground  surface 
and  the  atmosphere;  this  point  is  illustrated  by  the  fact  that  better  results 
can  be  obtained  by  using  17  channels  than  7  channels. 

The  last  four  columns  in  Table  2  show  results  in  the  retrieval  process 
as  the  group  of  coefficients  [Fp  are  smaller  by  10  percent  (represented  by 
F^q  in  the  table)  than  the  real  values  in  the  retrieval  process.  We  can  see 
that  there  are  very  small  differences  from  results  obtained  by  use  of  the 
correct  coefficient  [F^]  under  the  corresponding  conditions.  Only  when  both 
c  and  Ts  are  unknown  are  there  somewhat  greater  errors  (respectively  at  1.79*C 
and  0.0237)  for  the  obtained  retrieval  ground-surface  temperatures  and  thermal 
radiation  efficiencies;  hwoever,  the  effect  of  the  retrieval  atmospheric 
temperature  is  not  great  even  in  this  situation.  This  illustrates  that  the 
radiation  term  reflected  to  space  from  the  ground  surface  is  small  compared 
with  the  emission  term.  Therefore,  some  errors  will  not  cause  a  serioi  s 
effect  on  retrieval  in  the  calculation  coefficient  [F! ] .  This  point  is 
advantageous  in  actual  application. 

In  the  situation  when  the  atmospheric  pressure  at  the  ground  surface  is 
some  other  value,  similar  results  can  also  be  obtained.  Figure  2  shows  results 
when  the  ground- surface  atmospheric  pressure  Ps=700  mb.  We  can  see  that  quite 
good  results  can  still  be  obtained  when  both  the  ground-surface  parameter*, 
c  and  Ts,  are  unknown.  Figure  2  (a)  shows  an  example;  Figure  2  (b)  shows  the 
root  mean  square  deviations  of  all  layers  below  10  mb  as  obtained  in  289 
numerical  experiments. 

Since  the  situation  of  the  ground  surface  varies  with  locations  and 
seasons,  the  thermal  radiation  efficiency  at  the  ground  surface  varies  with 
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[Fig.  1  continued] 

Key:  (a)  Pressure  at  particular  layer;  (b)  Ground  surface;  (c)  Temperature; 
(d)  Initial  value  of  temperature  contour  line;  (e)  Retrieval  results  when 
both  e  and  Ts  are  unknown;  (f)  Retrieval  results  when  both  e  and  Ts  are  known; 
(g)  "Real"  temperature  contour  line;  (h)  Mode  of  atmosphere;  (i)  Guangzhou. 
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Fig.  2.  Temperature  retrieval  results  when  the  ground- 
surface  atmospheric  pressure  Pss700  mb. 

Key:  (1)  Pressure  (mb)  at  the  particular  layer;  (2) 
Ground  surface;  (3)  Temperature;  (4)  Initial  values  of 
temperature  contour  lines;  (5)  Retrieval  results  when 
both  c  and  Ts  are  unknown;  (6)  "Real"  temperature 
contour  line  (e=0.8S)  with  ground-surface  atmospheric 
pressure  Pss700  mb;  (7)  Solid  line  shows  results  when 
[Key  continued  on  following  page] 
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[Key  (of  Fig.  2)  continued] 

both  £  and  Ts  are  unknown;  the  root  mean  square 
deviation  of  all  layers  is  2.19*C;  (8)  Circles  (*) 
indicate  results  when  both  e  and  Ts  are  known;  the 
root  mean  square  deviation  of  all  layers  is  2.10*C. 


wavelengths.  Therefore,  it  is  a  very  difficult  problem  to  accurately  determine 
(in  advance)  the  ground -surface  parameters,  e  and  T$,  by  some  means  while 
seeking  retrievals  of  the  atmospheric  temperature.  However,  from  the  above- 
mentioned  results,  the  ground  surface  can  be  considered  as  a  gray  body  in  a 
spectral  region  of  the  selected  group  of  temperature-measurement  channels 
(actually,  within  the  range  of  several  channels  when  the  ground- surface  trans- 
missibility  is  not  equal  to  zero).  Thus,  no  other  means  is  required  to  know, 
in  advance,  the  ground-surface  parameters;  the  atmospheric  temperatures  and 
the  ground- surface  parameters  can  be  simultaneously  obtained  directly  fron 
these  temperature-measurement  channels.  This  is  very  significant. 

The  author  expresses  his  gratitude  to  Comrade  Zhou  Xiuji  for  his  revision 
and  for  providing  valuable  counsel. 
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Presonted  In  this  paper  Is  a  formula  of  oaloulating  brightness  temperature  in  short* 
wave  IR  atmospheric  window  from  instantaneous  radianoe  measurements  of  two 
channels  in  3.7~4.t>^m  spectral  band  by  polynomial  tit  technique.  Its  accuracy  is 
within  ±0.4  K.  The  approach  for  determining  oorreotions  of  retieoted  solar  radianoe  on 
measurements  in  4.8jim  OO,  absorption  band  is  briefly  described. 


I.  FOREWORD 

There  are  two  channels  (HIRS/2-18  and  HIRS/2-19)  of  a  short  wave  infrared 

(IR)  atmospheric  spectral  window  in  the  improved  model  of  high  discriminability 

IR  detector  (HIRS/2)  in  the  TIROS-N  meteorological  satellites  of  the  United 

States;  the  number  of  center  waves  is,  respectively,  25 11. 95  cm  1  and  2671.18 

cm"1,  marked  as  W.  and  W,.  8y  using  the  simultaneous  radiance  measurements  of 
1  1  28 


these  two  channels,  the  reflectivity  of  the  earth's  surface  can  be  calculated, 
and  the  effect  on  the  temperature  detection  channel  of  the  4.3  ym  CO.,  absorption 
band  on  the  solar  radiation  v'eflected  from  the  earth's  surface  can  be  estimated. 

II.  CALCULATION  KETHODS 


The  response  to  water  vapor  by  the  short  wave  IR  atmospheric  spectral 
window  is  considerably  weaker  than  that  of  the  long  wave  IR  atmospheric  spectral 
window  (11.0  ua) .  Therefore,  not  as  in  the  case  of  the  long  wave,  water  vapor 
absorption  is  a  factor  that  should  be  considered  first.  However,  the  extra¬ 
terrestrial  radiation  reflected  from  the  earth's  surface  should  be  considered; 
the  reflective  component  of  the  downward  radiation  to  the  atmosphere  (in  the 
reflected  radiation)  is  less  than  1  percent  of  the  radiation  at  the  earth's 
surface;  therefore,  the  reflective  component  can  be  neglected.  However,  the 
reflected  solar  radiation  is  about  1/7;  this  is  the  factor  having  a  major 

.HI 


effect  on  short  wave  IR  remote  sensing  during  daytime1 


For  the  short  wave  IR  atmospheric  spectral  window,  the  radiation  transmis¬ 
sion  equation  can  be  simplified  in  [2]  under  cloudless  atmospheric  conditions: 


Jt(H’)- r,/}(fr,  T,)+rj0m^ 


(1) 


In  the  equation,  R(W)  is  the  measured  radiation  intensity  when  the  intermediate 
and  short  wave  IR  atmospheric  spectral  window  channel  of  the  IR  detector  (the 
number  of  center  waves  is  indicated  !•/  W)  aims  at  the  globe.  The  first  term  on 
the  right  side  of  equation  (1)  represents  the  contribution  to  R(W)  by  solar 


radiation  reflected  from  the  earth's  surface.  ys  an*l  es  are*  respectively. 


the  reflectivity  and  thermal  radiation  efficiency  of  the  earth's  surface; 


is  the  cosine  of  the  zenith  angle  e0  of  the  sun;  when  gQ  is  greater  than  90*, 


uQ*0.  B(v,  T)  is  the  spectral  radiation  intensity  (also  called  the  Planck 


function)  at  wave  number  v  of  a  black  body  with  temperature  T. 


B(r,  7,'r- 


«V* 


exp  U »r 


~T)  - 1  ’ 


In  the  equation,  Cj  and  represent,  respectively,  the  first  and  second 
radiation  constants;  c^hc4-  and  c,»hc/k. 


I  is  the  solar  radiation  intensity  arriving  at  the  earth's  surface; 
o 

the  approximate  expression  of  is: 

7#(IK)  -  B(W,  (2a) 

In  the  equation,  Tgun  is  the  brightness  temperature  of  the  sun,  generally  at 

580(H6000K;  n  is  the  included  solid  angle  of  the  earth  by  the  sun  and  is  a 
sun  . 

constant,  approximately  7x1(7  sr.  Here  we  assume  that  in  magnitude  the 
direct  solar  radiation  is  equivalent  to  radiation  intensity  of  a  black  body 
with  temperature  T 

sun 

If  the  microabsorption  effect  of  the  gaseous  components  of  the  atmosphere 
on  the  solar  radiation  in  the  short  wave  IR  spectral  window  is  included, 
equation  (2a)  is  revised  as: 

/.(JK)-Il(W'.  *■♦■*"( WO.  (2b) 


In  the  equation,  t (W)  is  the  transmissibility  at  the  upward  perpendicular 
direction  of  the  local  site;  eQ  is  the  zenith  angle  of  the  sun;  and  e  is  the 
zenith  angle  of  the  satellite.  In  these  expressions,  the  form  of  equation  (I) 
remains  unchanged. 

By  using  Smith's  assumption ^  that  Ys  and  es  are  not  related  to  the 
approximate  value  of  the  wavelength  in  a  narrow  short  wave  IR  spectral  region 
from  3.7  um  to  4.2  ym,  we  obtain  equation  (3)  from  equations  (1)  and  (2): 

T.)  __  By  MY  T.J\ 

B^Wl)-t.BKWl,  T,)  ~B(WX,  J'm)  9  w 

For  HIRS/2-18  and  HIRS/2-19,  if  Tsun»5800K,  the  right  side  of  equation  (3) 
becomes: 


B(  W,, 

WWi.  T~) 


£'-1.10035, 


(4) 
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According  to  the  definition,  the  brightness  temperature  Tg  satisfies 
the  following  equation: 

•.BW,  (5) 

Equation  (3)  can  be  rewritten  as: 

T,)- K D(WU  Tm)  -  -  k'R(Wt)0  (6) 


fl 


Equation  (6)  is  a  transcendental  equation  related  to  Tg.  We  used  the 
polynomial  fit  method  to  establish  the  apparent  relationship  between  Rfwp 
on  one  hand,  and  Tg  and  the  radiation  measurement  value  R(Wj)  on  the  other. 
Let 

/-JKVF^-ATJKH',).  (7) 


We  take  the  cubic  polynomial;  i.e., 

r.-ga.tln(-/)3‘,  W 

From  equations  (6)  and  (7),  we  can  calculate  the  coefficient  a^  first.  By 
selecting  the  temperature  range  from  200K  to  340K  (corresponding  to  two 
extreme  temperatures  of  the  earth),  the  least  squares  method  is  used  (on  two 
short  wave  1R  atmospheric  spectral  window  channels  of  TIROS-N)  to  establish 
a  linear  equation  set  concerning  a^  and  to  adopt  a  main  element  cancellation 
method  for  solving  the  equation  set.  Table  1  shows  the  calculation  results. 

Table  1.  Coefficients  of  the  calculated  brightness  temperatures  of  HIRS/2 
short  wave  IR  spectral  window  in  TIROS-N. 
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Table  2.  Calculated  value  (T.)  and  simulated  value  (T)  by  use  of  polynomial 
fit  method.  b 


210 

240 

280 

210.08 

228.95 

259.88 

279.97 

0.08 

-0,05 

-  0.14 

-O.IO 

It  should  be  pointed  out  that  although  the  effect  of  water  vapor  and  other 
gaseous  components  is  relatively  small  in  the  short  wave  IR  atmospheric  spectral 
window,  their  influence  on  the  measurement  values  should  still  be  considered, 
especially  in  the  case  when  the  water  vapor  content  in  the  low  atmospheric 
layer  is  relatively  great  (in  this  paper,  it  is  assumed  that  corrective 
measures  are  made  on  absorption  of  these  components  before  application  of  equa¬ 
tion  (8)). 

III.  CORRECTIONS  OF  REFLECTIVITY  (TERM)  OF  THE  EARTH’S  SURFACE  AND  THE 
REFLECTED  SOUR  RADIATION  TERM 

After  calculating  the  brightness  temperature  Tg  by  using  the  measurement 
values  of  IR  atmospheric  spectral  window  channels  by  use  of  the  satellite-borne 
IR  detector,  the  reflectivity  ys  is  (based  on  equation  (1)): 

v  _  T ,) 

7t  ' 

If  W  is  substituted  by  Wj,  we  obtain  by  using  equation  (2b): 


„  _  r.)  ro~..  i 

r* - HWT7Z)  1  (  J  • 


The  corrected  term  SC(V)  of  the  reflected  solar  radiation  of  a  certain 
channel  (indicated  by  v,  the  number  of  center  waves)  within  the  CO,  absorption 
zone  of  4.3  urn,  actually  is  the  second  term  on  the  right  side  of  equation  (1). 


SC(«0 


kop.)- W..  r„) 

r  r(«)  V"*4-* 

b\Wu  T„) 

L  t^iY-1 

(10) 


In  the  equation,  the  zenith  angle  8o  of  the  sun  can  be  calculated  by  using  the 
astronomical  formula  at  the  pltce  of  the  satellite's  observation  time.  There 
is  another  parameter,  the  transmissibility,  which  should  be  derived  by  using 
an  appropriate  method  in  advance. 
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SEVERAL  PROBLEMS  IN  DEVELOPING  SATELLITE  BORNE 
INFRARED  REMOTE  SENSOR  AND  WAYS  TO  SOLVE  THEM 

Ki  ng  IIuixing 

(Shanghai  hutUute  Of  Technical  Phytic t,  Academia  Stiuca) 


ABSTRACT 

Several  spoolally-considoied  problems  in  developing  a  satellite  borne  infrared  re¬ 
mote  sensor  are  discussed.  These  problems  include:  tbo  definition  of  scanning  period, 
the  soanner  and  its  lubrication  in  outer-spaoo  environment  of  very  high  vaouum,  and 
infrared  radiation  calibration  and  so  on.  The  ways  to  solve  thorn  are  suggested. 
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INFRARED  REMOTE  SENSING  Cf  OVER- LAND 
ATMOSPHERIC  TEMPERATURE  PROFILES  FROM  SATELUTE 

Zhao  Gaoxiano 

(Institute  o]  Atmo'phmc  thysics,  Acmiemis  Sir.ua) 


ABSTRACT. 

The  effects  of  the  surfaoo  p.ir?  motors,  t  and  T ,,  on  alraosphorio  temperature  rotrie- 
rale,  and  the  possibility  of  recovering  simultaneously  both  t,  T,  and  otmoepherio 
temperature  profiles  are  studied.  The  results  of  numorioal  simulations  show  that,  when 
me  of  the  two  surfaoo  parameters  is  known,  atmospherio  temperature  profiles  oan  be 
lerived  and  at  the  sumo  titno  the  other  surfaoo  parameter  can  be  determined  aoourately; 
ind  even  if  both  t  and  T,  are  unknown,  atmosphorio  temporatuxe  profiles  and  the  two 
parameters  oan  still  be  recovered  with  good  aocuraoios. 
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A  METHOD  OF  DETERMINING  THE  BRIGHTNESS 
TEMPERATURE  IN  SHORT-WAVE  IR  ATMOSPHERIC  SPECTRAL 
WINDOW  FROM  SATELLITE  MEASUREMENT 


Wang  Qinmo 

{Satellitt  Uettoroloj/icai  Ctnler.  CU.B) 

ABii’RACT 

Presented  in  this  japor  is  a  formula  of  calculating  brightness  temperature  In  short* 
wave  lit  atmospheric  window  from  instantaneous  radianoe  measurements  of  two 
channels  in  3.7~4.l)^m  spootral  band  by  polynomial  lit  toohnique.  Its  aoouraoy  is 
within  ±0.4  K.  The  apprcaoh  for  determining  corrections  of  reHeoted  solar  radianoe  on 
measurements  in  4.?,/x m  CO,  absorption  band  is  briefly  described. 
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